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Surface properties of the clean and Au/Pd covered Fe304(lll): 

a DFT and BFT+U study 

Adam Kiejna0 Tomasz Ossowski, and Tomasz Pabisiak 
Institute of Experimental Physics, University of Wroclaw, Plac M. Borna 9, PL-50-204 Wroclaw, Poland 

The spin-density functional theory (DFT) and 'DFT+U with Hubbard U term accounting for on- 
site Coulomb interactions were applied to investigate structure, stability, and electronic properties of 
different terminations of the Fe304(lll) surface. All terminations of the ferrimagnetic Fe304(lll) 
surface exhibit very large (up to 90%) relaxations of the first four interlayer distances, decreasing 
with the oxide layer depth. Our calculations predict the iron terminated surface to be most stable 
in a wide range of the accessible values of the oxygen chemical potential. The adsorption of An and 
Pd on two stable Fe- and O-terminated surfaces is studied. Our results show that Pd binds stronger 
than An both to the Fe- and O-terminated surface. DFT+U gives stronger bonding than DFT. The 
bonding of both adsorbates to the O-terminated magnetite surface is by 1.5-2.5 eV stronger than 
to the Fe-terminated surface. 



PACS numbers: 68.35.-p, 68. 43. Be, 68.47.Gh, 73.30.+y 



I. INTRODUCTION 

Magnetite (Fe304), the earliest known magnetic mate- 
rial, is an abundant mineral in nature and is important 
due to its potential applications in magnetoelectronics 
and in catalysis. It is one of the most interesting iron 
oxides formed during the corrosion processes (rusting). 
Ultra-thin films and nanostructures formed by noble met- 
als on iron-oxide support show enhanced catalytic prop- 
erties compared with a clean oxide surfaces. Thus, under- 
standing the properties of magnetite surfaces is of utmost 
importance from the viewpoint of basic science and ap- 
plications. However, even clean iron oxide surfaces are 
relatively little explored which is connected with difficul- 
ties in the preparation of well defined surfaces [l|, |3| • 

At room temperatures and under normal pressure con- 
ditions magnetite crystallizes in the inverse spinel struc- 
ture (space group Fd3m) in which tetrahedral positions 
are occupied by ferric (Fe'^+) while octahedral ones con- 
tain equal number of ferric and ferrous (Fe^+) iron atoms. 
Magnetite's primitive rhombohedral cell contains two for- 
mula units of Fe304 and its volume is equal to one quar- 
ter of the spinel unit cell. Bulk magnetite is a semi-metal. 
Since the Fe^+ ions are aligned antiferromagnetically and 
the ratio of Fe'^^ and Fe^+ ions is 2:1 the overall crystal 
structure is ferrimagnetic. At T « 120 K, the metal- 
insulator Verwey transition occurs which is connected 
with a long-range change of the degree of localization 
of electrons in the octahedral Fe atoms jl-dl ■ 

The structure of magnetite can be also represented by 
the hexagonal conventional unit cell which contains eight 
formula units. In this stacking, (111) oriented layers of 
oxygen atoms separate alternate Fe monolayers of octa- 
hedral (Fcocti) sites and Fe trilayers consisting of a Feoct2 
monolayer in octahedral sites with Feteti and Fetet2 layers 
in tetrahedral sites on either side. Thus the stacking se- 
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quence of the atomic planes perpendicular to the [111] 
direction can be written [l| as Feoct2-Feteti-0i-Feocti- 
02-Fetot2-- The (111) surface, which is the dominant 
cleavage plane of magnetite and is often exposed on nat- 
urally grown crystals can have six different terminations 
[l| . Only four of them have been confirmed experimen- 
tally, namely, Feoct2 (3], Fcteti 13, and densely packed 
oxygen planes O2, and Oi d, U Both iron and oxy- 
gen terminated Fe304(lll) surfaces are polar of type III 
according to Tasker's classification [l^]- The surface ter- 
mination is very sensitive to the preparation conditions 
of the samples [ij which are either Fe304 single crys- 
tals or epitaxial Fe304 films on single crystal substrates. 
Despite many experimental investigations, information 
about the (111) surface at the atomic level is very scarce. 
There are not many data about its magnetic properties 
or how its structure and composition change with tem- 
perature and oxygen pressure. Magnetite is a strongly 
correlated system and DFT with standard (local or semi- 
local) exchange-correlation functionals does not allow for 
a correct description of its electronic structure because of 
inadequate treatment of the strong Coulomb interaction 
between 3d electrons localized on the Fe ions. This short- 
coming of DFT is corrected in practice by either of two 
semi-empirical approaches: hybrid functionals, where the 
exact Hartree-Fock exchange is partially mixed with the 
DFT exchange, or DFT-I-J7, where the on-site Coulomb 
repulsion is described by an additional Hubbard term U. 

The question what is the most stable termination of the 
Fe304(lll) surface is stifi open and controversies about 
its structural details seem to remain. An earlier scanning 
tunneling microscopy (STM) study [61] reported two co- 
existing surface terminations. One was assigned to Fcocti 
atoms and the other to Feoct2~Fetoti layers. A low energy 
electron diffraction (LEED) study [21 concluded that the 
(111) surface termination corresponds to 1/4 monolayer 
of Fe atoms over a hexagonal close-packed O layer under- 
neath. This disagreed with ab initio periodic Hartree- 
Fock calculations of Ahdjoudj et al. 11] who found the 
Feoct2~Fetoti bilayer to be the most favorable termina- 
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tion of the clean surface. Lemire et al. [l2| studied the 
surface structure of Fe304(lll) films by CO adsorption 
and concluded that the (111) surface is terminated with 
Feoct2- Recent full-potential linearized augmented plane- 
wave (FP-LAPW) calculations by Zhu et al. pJi] deter- 
mined the structure, composition and relative stability of 
five different terminations of the Fe304(lll) surface. The 
effect of different functionals for exchange and correlation 
energy (generalized gradient approximation (GGA), local 
density approximation-|-[/ (LDA-I-C/)) has been also dis- 
cussed. According to that study the Feoct2 termination 
is the most stable one. This, however, was not confirmed 
by a more recent first principles calculations performed 
by Grillo et al. [11] within GGA+U and Martin et al. 
[l5| using GGA, according to which the Fctcti-terminated 
surface has the lowest surface energy. Recent STM exper- 
iments have reported that the stoichiometric (111) sur- 
face of magnetite corresponds to Fctcti termination p^ . 
Another, recent combined STM and first principles cal- 
culation study [l3l also predicted the Fctoti termination 
as the most stable one. 

The adsorption of metal atoms on magnetite was very 
rarely studied. On the theoretical side only the inter- 
action of alkali metal atoms with the Fctoti-terminated 
Fe304(lll) surface has been studied using DFT at the 
GGA level [18]. To our knowledge, studies of the adsorp- 
tion of noble and transition metal atoms on Fe304(lll), 
which are important in catalysis, have not been reported 
so far. 

In this work we revisit first the calculations for different 
terminations of the clean Fe304(lll) surface using DFT 
and DFT-I-C/ approaches, in order to explore the infiuence 
of strong on-site electronic correlations on the structure 
and physical properties of magnetite, and to form a firm 
basis for our studies of Au and Pd atom adsorption on 
the Fe304(lll) surface, which is the main subject of this 
work. 



II. METHODS AND COMPUTATIONAL 
DETAILS 

The calculations presented in this work are based on 
the spin density functional theory as implemented in the 
VASP package 19, 20]. The calculations employed the 
GGA-PW91 version of the exchange and correlation en- 
ergy functional |21)] with the spin interpolation of Vosko 
et al. HI]) and the GGA plus on-site Coulomb interac- 
tion term U (GGA-I-J7) using the Dudarev et al. [l^ 
approach. Following previous calculations [1, HI] for the 
bulk and surfaces of Fe304 the GGA-I-J7 calculations 
were performed with the effective parameter of interac- 
tion between electrons UcS = U — J = 3.61 eV (the 
Coulomb and screened exchange parameters {U, J) ~ 
(4.5,0.89) eV, respectively). The electron ion-core inter- 
actions were described by the projector augmented wave 
(PAW) method [1^. The plane waves basis with cut-off 
energy of 500 eV and the conjugate gradient algorithm 




FIG. 1: (Color online) The Fe304(lll) slab used in the sur- 
face calculations. The successive terminations were created 
by removing top and bottom layer from the thickest Fetoti 
terminated slab. Iron and oxygen atoms are represented by 
small and large balls respectively. The right hand side figures 
show top views of the considered terminations. The parallel- 
ograms mark the 1x1 surface cell applied in the calculations. 



were applied to determine the electronic ground state. 
The integrations over the Brillouin zone were performed 
using Monkhorst-Pack grids [1^ . A Gaussian broadening 
of the Fermi surface of 0.2 eV was applied to improve the 
convergence of the solutions. The results presented in 
this work were obtained using fc-point meshes of 6x6x6 
for the bulk, and 6x6x1 for surface calculations, which 
allowed to obtain total energy convergence within 1 meV. 
For surface calculations F-centered grids were used. 

Calculations of the Fe304 bulk structure showed that 
among the nonmagnetic and different magnetic phases, 
the ferrimagnetic phase is most stable, with the magnetic 
moments on the Feoct atoms antiparallel to those on the 
Fctct atoms. The lattice constant and bulk modulus of 
the magnetite crystal calculated within GGA, 8.377 A 
and 172 GPa, respectively are in ver y g ood agreement 
with experiment (8.396 A |23K 8.393 A [Ig; 181 GPa 0) 
and other GGA calculations |l5l [29j . The GGA magnetic 
moments on the Fctct and Feoct atoms are — 3.45/iB and 
3.49-3.61/iB respectively. The moments on the O atoms 
are much smaller (^0.08/is). The total magnetic mo- 
ment (3.66/iB) per formula unit, is about 0.4/iB lower 
than the experimental value [s^]- The GGK+U calcula- 
tions give the lattice constant (8.473 A) and bulk mod- 
ulus (182 GPa) in good agreement with experimental 
data and other calculations [l3, [13 ■ They improve the 
magnetic moments on the Fe atoms compared to GGA 
and give 4.04^^, and 3.91-3.95^b on the Fctot and Fcoct 
atoms, respectively. With the moments on the O atoms 
reduced by 50%, the total magnetic moment per formula 
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FIG. 2: (Color online) The dependence of surface energy on the oxygen chemical potential, fio, showing relative stability of six 
terminations of the Fe304(lll) surface calculated within the GGA and GGA+(7. Vertical, dashed lines mark allowed range of 
the oxygen chemical potential. 



unit is 3.97ub, in very good agreement with experiment 

(4.05/zs m)- 

The optimized lattice parameters were used to con- 
struct the Fe304(lll) surface slabs consisting of 19 to 
29 atomic layers separated by a vacuum region. Start- 
ing from the thickest symmetric slab for the Fetcti- 
termination (Fig. [T|) the other terminations were created 
by striping off subsequent atomic layers, without chang- 
ing the supercell size, both from the top and the bottom 
of the slab. Thus the surface was separated from its pe- 
riodic replicas by vacuum region ranging from 15.5 to 24 
A (the latter is for the thinnest Feoct2-terminated slab). 
In surface calculations the positions of all atoms were 
relaxed until the forces were smaller than 0.02 eV/A. 

The stability of different terminations of the 
Fe304(lll) surface as a function of oxygen partial pres- 
sure was considered based on the ab initio thermody- 
namics jm in which the surface free energy ^{T,P) is 
expressed as a function of the pressure and temperature 
by the chemical potentials /ipo, A'o of the constituents: 



7(r,P)^^[Gito. 



7VfoMfc(T,P) -iVoMo(T,P)] , 

(1) 

where the Gibbs free energy, Gp^^'j-,^, can be expressed 
by the total energy of the slab Ef^^, and iVpo and iVo 
represent the number of Fe and O atoms in the system. 
The Gibbs free energy per magnetite formula unit, g, 
is related to the chemical potentials of iron and oxygen 
through the relation ^p^^'q^ = 3/iFc-l-4^0- Consequently, 
the surface energy as a function of oxygen chemical po- 
tential can be written as 
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(2) 

where the Gibbs free energy is approximated by the inter- 
nal energy from DFT calculations [Slj. /xq is referenced 
with respect to the chemical potential of oxygen in a gas 
phase = ^£"02'' '^h^i'^ the total energy of an oxygen 



molecule i?Q°* is calculated in a large box. 

Au and Pd atoms were adsorbed in four different ad- 
sorption sites of the 1x1 surface unit cell, on both sides of 
the relaxed, symmetric Fe304(lll) slab. The adsorption 
binding energy was calculated as 

£;^d = -(£;X^"b_^sub_2£;X)/2, (3) 

where £:X/sub jg ^j^g ^^^^^^ energy of the slab covered with 
adsorbate X, E^^'^ represents the energy of the relaxed 
bare oxide support, and is the energy of a free adsor- 
bate atom. 



III. RESULTS 

A. Clean Fe304(lll) surface 

The variation of surface energy as a function of the 
oxygen chemical potential /iq is displayed in Fig. [2j The 
accessible range of /io is limited, at the lower limit by 
the onset of magnetite decomposition to form bulk iron 
and, at the upper limit, by start of oxygen condensation 
on the magnetite surface. It is seen that surface compo- 
sition of the Fetoti termination minimizes the surface free 
energy over a wide range of for both types of calcula- 
tions, in agreement with previous predictions [l3. fl6l [TtII . 
However, at low which correspond to very low oxy- 
gen pressures the GGA-|-?7 results show that the Feoct2 
termination may turn out to be stable. Zhu et al. [l3l | 
have predicted this termination to be most stable. The 
GGK+U does not give the stable oxygen terminations 
of the Fe304(lll) surface seen at higher pressures in the 
GGA results. In the following we investigated only four 
clean surfaces of magnetite - the two iron terminated 
(Fetoti, Feoct2), and two oxygen terminated (Oi and O2) 
ones. 

After structural optimization the surfaces are strongly 
relaxed. The relaxation of the interplanar distance is cal- 
culated as Aij = {dij — d)/d, where is the distance 



4 




e 



u 

Pi 30 

-30 
-60 



n I ^ \ \ \ \ I I \ 




~i — r 
Fe 



oct2 



/fe 



J L 



FIG. 3: (Color online) Relaxations, Aij, of the interplanar 
distance for different Fe304(lll) terminations. Experimental 
data for the Feteti termination are taken from Ref. 0. In the 
ideal crystal bulk the subsequent separations of the (111) lay- 
ers calculated within GGA (GGA+f/) are as follows: Fetoti- 
Oi, 0.637 (0.644) A; Oi-Feocti and Feocti-02, 1.176 (1.190) 
A; 02-Fetet2, 0.637 (0.644) A; Fe tet2-Feoct2 and Feoct2-Fetoti , 
0.605 (0.611) A. 



of the relaxed i and j — i + 1 planes, and d is the cor- 
responding distance in the bulk. The calculated relax- 
ations (Fig. [3]) show an oscillatory character, though not 
very regular, of the contraction-expansion type. In gen- 
eral, the geometry changes calculated within the GGA 
and GGA-I-J7 approaches are qualitatively similar and 
predict a large contraction of the first interplanar dis- 
tance. They differ, however, in the topmost layer re- 
laxation of the Feoct2 termination, where GGA predicts 
a large (40%) contraction of the Feoct2-Fetoti-layer dis- 
tance, whereas GGA+U shows no relaxation. Top layers 
relaxation calculated within LDA+U are substan- 
tially larger compared with our GGA values. As seen 
from Fig. [3] for the Fctoti terminated surface the relax- 
ations are in very good agreement with experimental data 
0, and agree well with results of the LDA+J7 calcula- 
tions fis!]. The largest relaxations (up to about 90%) are 
observed in the 02-termination. The first and second in- 
terlayer distances exhibit large contractions of about 50% 
GGA (53%, GGA+U) and 56% GGA (60%, GGA+U), 
respectively. These are almost compensated by a large 
expansion of the third interlayer spacing (92%, GGA; 
85%, GGA+U). The relaxations of the O2 surface are 



much larger (up to three times) than those found from 
LDA+U They lead to a specific final configuration: 
a surface region of reduced thickness formed by a triple- 
layer, slightly separated from the deeper lying oxide lay- 
ers. Similar structures are observed for the other consid- 
ered terminations: separated triple layers for the Fcteti- 
and Feoct2- terminations, and a separated double layer for 
the Oi-terminated surface. In, general, the relaxations of 
the 0-terminated surfaces of magnetite are much larger 
than those observed on 0-terminated hematite (0001) 
surface [32]. The topmost layer relaxations on the iron 
terminated surface of magnetite (111) show similarity to 
those calculated for hematite (0001) surfaces [13 • The 
calculated contraction of the spacing between the top- 
most Fe and O surface layers might be considered as a 
mechanism which reduces the total electrostatic dipole 
moment and thus stabilizes these polar surfaces. The 
ideal bulk 0-layers show only a small corrugation. On 
the Oi plane one of the O atoms is slightly shifted down 
giving a 0.03 A corrugation with respect to the average 
position of the plane. A similar corrugation is caused by 
an upward shift of one of the O atoms of the O2 plane. 

Figure m shows the calculated density of states (DOS) 
of bulk magnetite and the local density of states (LDOS) 
at the four Fe304(lll) surface terminations. As is seen 
the DOS and LDOS resulting from the respective GGA 
and GGA+U calculations differ substantially. Both 
methods predict half-metallic behavior [s^ of the bulk 
magnetite with the band gap in majority spin electrons 
which is much wider when calculated within GGA+U. 
However, all considered (but O2) surface terminations 
show metallic LDOS when calculated within GGA. The 
metallic character of the Oi termination agrees with the 
results of a combined spin-polarized STM and DFT study 
reported by Berdunov et al. Q. This prediction is 
at variance with the results from photoemission spec- 
troscopy studies [H, HI] and with our GGA-|-f7 results 
(Fig. 2]) which show that almost all terminations remain 
half-metallic. The exception is Feteti surface where the 
narrow (~0.3 eV) energy gap opens below the Fermi level. 
The location of the band gap alters from the majority 
to minority spin band, depending on the termination. 
Compared with the GGA the LDOS calculated within 
GGA+U extend over a wider energy range and the main 
weights of the Fctcti and Feoct2 states are shifted to lower 
energies. At the Feoct2 termination responsible for the 
half-metallic character are Fe Zd states, whereas at the 
Oi and O2 terminations half-metallicity is due to the hy- 
bridized O 2p and Fe 3d states. 

The difference in the electronic structure of the dif- 
ferently terminated surfaces is manifested in their work 
function values. It was calculated as the difference be- 
tween the electrostatic potential in the vacuum and the 
Fermi energy of the slab. The work function of the 
0-terminated surfaces is 2-3.5 eV (GGA) and 2-4 eV 
[GGA+U) larger than that of the Fe-terminated surfaces. 
The GGA work functions are about 0.4-1.0 eV lower 
than those determined from GGA+U (Table Gener- 
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FIG. 4: (Color online) Density of states of bulk magnetite and LDOS on atoms of six topmost atomic layers of different surface 
terminations decomposed into the contributions from the iron and oxygen atomic layers. Majority and minority spin states 
are respectively displayed as positive and negative. The left and right hand side panels show respectively GGA and GGA+t/ 
results. 



TABLE I: Work function of different Fe304(lll) surfaces. 



Termination 




Work function (eV) 


GGA 


GGA+f/ 


Feteti 


4.93 


5.48 


Feoct2 


3.58 


3.90 


Oi 


7.04 


8.09 


O2 


6.88 


7.66 



ally, the work function is lowest for the Feoct2-terminated 
surface and highest for the Oi termination. The calcu- 



lated Bader charges [36|, |37| on atoms at the exposed sur- 
face terminations show that the atoms in the topmost Fe 
layer, both of the Feteti and Feoct2-terminated surfaces, 
gain electron charge compared to that on atoms in the 
bulk crystal layers (0.28e; 0.76e and 0.27e, respectively). 
The oxygen atoms of the 0-terminated surfaces, as well 
as the 0-atoms of the Feteti termination, lose electrons. 
In contrast, in the Feoct2 termination, not only the Feoet2 
and Feteti atoms gain electrons but also most of the O- 
atoms of the subsurface layers except for one of the atoms 
of the Oi-layer which loses 0.14e. 

The changes in the magnetic moments on the iron 
atoms due to the presence of the surface and its re- 





FIG. 5: (Color online) Schematic illustration of the initial Au 
and Pd adsorption places at the Feteti terminated surface in 
a top (left) and side (right) view. 



laxation are limited practically to the 2-3 outmost sur- 
face layers. The directions of the magnetic moments do 
not change and the surface remains ferrimagnetic, but 
the magnitudes of the moments are smaller by about 
0.1-1.6/LtB (GGA) and 0.05-0. 56/LtB (GGA-f [/) compared 
to the corresponding atoms in magnetite bulk. The 
strongest change of the magnetic moments is observed 
for the topmost surface layer. At the Feteti-termination 
the moments of the topmost Fe layer are reduced to 
3.11/iB (GGA) and 3.51^3 {GGA+U). At the Feoct2- 
terminated surface the moment is enhanced to A.02fiB 
(GGA) and 4.45/xb (GGA+U). As can be seen from 
the lowest panel of Fig. |31 this is connected with the 
big change in the relaxation of the topmost Feoct2 layer 
which for the GGA+U is practically reduced to zero. The 
atoms of subsurface 0-planes of the Fe-terminated sur- 
faces exhibit small magnetic moment ^O.l-OAfiB- The 
GGA+U moments on oxygen atoms of the Oi and O2 
terminated surfaces are in the range 0.15-0.23/iB and 
(— 0.02)-0.20/Lt_B, for the Oi and O2 termination, respec- 
tively. For the O2 termination one of the O atoms in the 
topmost layer has a very small (~0.02/iB) local moment 
opposite to the other three atoms. The moments on the 
Fe atoms of subsurface planes in these two oxygen ter- 
minations show a significant reduction when calculated 
within GGA. In contrast, the GGA+U values deviate 
relatively little from the moments on bulk atoms. 



B. Adsorption on the Fetcti terminated surface 

The adsorbate atoms were initially placed well above 
four different sites (Fig. [5]) of the Feteti -terminated sur- 
face: in site A, on-top of the Fctcti atom, site D above the 
Oi atom, or in one of the hollow sites, B or C, formed re- 
spectively in a deep hollow over the O2 atom, and in the 
hollow above the Fcocti atom. All of the adsorption sites 
appeared to be stable after structure relaxations. The 
most stable position for Au adsorption is site A, whereas 
Pd atoms adsorb preferably in site C. 

The changes in the surface relaxation due to Au/Pd 
adsorption are presented in Fig. [B] The adsorption 
of Au and Pd in any considered adsorption site tends 
to suppress the surface relaxation of the first interpla- 
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FIG. 6: (Color online) Relaxation of the interplanar distance 
after Au/Pd adsorption in most stable sites on the Fetcti ter- 
minated Fe304(lll) surface. 



nar distance of the Fetcti termination compared to the 
clean surface. The relaxation patterns resulting from 
GGA and GGA+U calculations are very similar, with 
the relaxation of the first four layers being larger for 
GGA+U . For Au adsorbed in the most preferred on- 
top site A, a large contraction of the clean surface (35%, 
GGA; 40%, GGA+C/) is converted into distinct expan- 
sion (~10% GGA; ~15%, GGA+U), which makes the 
separation larger than respective one in the bulk crystal. 
A Pd adatom in the on-top site A causes much smaller 
changes of the interplanar distance which remains con- 
tracted compared to that in bulk magnetite. In general, 
the changes in the surface geometry due to Pd adsorption 
are much smaller than those caused by Au. A Au/Pd 
atom adsorbed in any other considered site does not 
change qualitatively the relaxation pattern. For both ad- 
sorbates, the distances between deeper oxide layers show 
relatively small changes. 

Figure [7] displays the calculated adsorption energy and 
adatom-surface distance for a single Au/Pd atom adsorp- 
tion. As is seen, in general, the GGA+U binding is 0.2- 
0.5 eV stronger than that resulting from GGA. The Pd 
binding to the Feteti terminated surface is stronger than 
that of Au in all sites, except for site A. This agrees with 
the results for Au/Pd adsorption on the iron terminated 
(0001) surface of hematite [321 • The adsorption energies 
in different sites are less differentiated for Pd than for Au. 
For Au adatoms the most stable position is site A (on top 
of the Fetcti atom), whereas for Pd it is the threefold O- 
coordinated hollow C (Fig. [S]). The binding energy of 
the Au (Table [TT]) in the most preferred on-top site A is 
1.66 eV (GGA) and 1.98 eV {GGA+U). In this site, the 
Au-Fctcti bond is perpendicular to the magnetite surface 
with a bond length of 2.43 A (GGA and GGA+U). For 
Pd adsorption the most stable site is the threefold O- 
coordinated hollow site C with an adsorption energy of 
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FIG. 7; (Color online) Adsorption energy, i?ads, work func- 
tion change, A<1>, and adatom-surface distance for a single X 
(=Au,Pd) atom adsorbed in different sites of the Feteti ter- 
minated surface. 



TABLE II: Adsorption energy, Sads, adatom-surface distance, 
dx-Foteti 1 ^iid work function change, A$, for X (=Au,Pd) 
adatom in most stable sites on the Fetcti-terminated surface. 
For each quantity the left and right hand side columns display 
GGA and GGA+U values, respectively. 



X (site) 


E,ds (eV) 


rfx- 


-Fototl (^) 


A$ (eV) 


Au (A) 


1.66 1.98 


2.41 


2.43 


0.84 0.52 


Pd (C) 


1.76 2.20 


1.06 


1.19 


-0.67 -1.20 



1.76 eV (GGA) and 2.20 eV (GGA+U). The respective 
Pd-0 bond lengths are 2.18 A, 2.46 A, and 2.46 A, which 
is not much larger (2.68 A) than the distance between the 
Pd and Fcocti atom in the third surface layer. 

The stronger binding of Pd than Au can be understood 
by inspecting the results of the layer-resolved LDOS pre- 
sented in Fig. [SI They were calculated for the Au in 
site A and the Pd in site C, which are their most sta- 
ble positions. The GGA LDOS of surface layers with Au 
and Pd show metallic character. The principal peaks of 
the Pd 4d electron states are higher in the energy and 
hybridize with the O 2p states closer to the Fermi level 
than the Au 5d peaks, which results in stronger binding 
of the Pd. In the case of Au adsorption the metallicity of 
the oxide is enhanced by the Au 5d states in the energy 



range between -3 eV and the Fermi level and a smaller 
contribution from the Au 6s states at the energies close to 
the Fermi level. In the GGA-I-J7 approach the adsorbed 
Au introduces a small density of 6s states at the Fermi 
level which contributes to the reactivity of this surface 
and changes the semiconducting character of the surface. 
The Pd 4c? states tend to close the energy gap on this 
oxide termination and convert it from semiconductor to 
semi-metal. For the GGA+U the main weight of the Pd 
4d states is about 1 eV closer to the Fermi level than that 
of the Au 5c? states, which again means stronger binding 
of Pd than Au. 

The differences between Au and Pd adsorption are also 
seen in the work function changes (Fig. [7]). In all consid- 
ered sites Pd lowers the work function while Au distinctly 
increases it only when adsorbed in the A site. In the re- 
maining sites Au either slightly reduces the work function 
(GGA+U) or has only a very small effect on it (GGA). 
In general, GGA predicts a higher work function increase 
(by ^0.3 eV) than GGA+U for Au adsorption and a 
smaller (by 0.4-0.7 eV) lowering of the work function due 
to Pd adatom. For Au adsorption in the most preferred 
site A, the work function is increased by 0.84 eV (GGA) 
and 0.52 eV (GGA+U). In contrast, the work function of 
the Pd/Fctcti system, with Pd in the most stable site C, 
is lower by 0.67 eV (GGA) and 1.20 eV (GGA+U) com- 
pared with that of the clean surface. The above work 
function changes are consistent with the electron charge 
transfer to/from the surface atoms due to adsorption of 
Au/Pd. The Bader charge analysis [H, [s^ shows that a 
Au adatom in site A gains electrons (0.27e, GGA; 0.32e, 
GGA+U) at the expense of the surface Fe and O atoms. 
In the case of GGA+U most of the charge (0.22e) is do- 
nated by the Fctcti atom from beneath the Au. In the 
case of Pd adsorption in the site C, the charge (0.42e, 
GGA; 0.36e, GGA+U) is transferred from the adatom 
to the surface, mostly to the Fetoti atom (0.13e, GGA; 
0.05e, GGA+U). 

The changes of the magnetic moments caused by 
Au/Pd adsorption are small and limited to substrate 
atoms in the 2-3 topmost layers. For Au in the most sta- 
ble site A, the GGA+U magnetic moment on the topmost 
Fcteti atoms is enhanced by about 0.42/xb with respect 
to that of the clean surface. For Pd GGA-t-J7 adsorp- 
tion in the most preferred site G the change is negligible 
(fnO.OlfiB). Small magnetic moments appear on the Au 
and Pd adatoms. The moment on the Au atom in dif- 
ferent sites is in the range of 0.14-0.20/iB (GGA) and 
0.10-0.19/,iB (GGA+U), the lowest values corresponding 
to the Au atom in most stable site A, and has the same 
direction as the moments on the Fe atoms in the topmost 
Feteti layer. The magnetic moment on the Pd adatom in 
sites A, C, and D is in the range of 0.07-0.27/i_b (GGA), 
and 0.03-0. 16/.tB (GGA+U), and has the same direction 
as the moment on the Fe atoms in the topmost surface 
layer. The magnetic moment on Pd at the site B is 
0.41^5 and 0.62^^ for GGA and GGA-I-J7, respectively, 
and is oriented upwards. The moment on the Pd adatom 
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FIG. 8: (Color online) Local density of states for Au and Pd adsorption in most stable sites on the Fetcti termination resulting 
from GGA and GGA+C/ calculations. LDOS of the adsorbate atom and of the three (two Fe and one O) topmost atomic layers 
of the magnetite (111) surface are shown. Corresponding LDOS for the clean Fetcti -termination are displayed for comparison 
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FIG. 9: (Color online) Schematic illustration of the initial Au 
and Pd adsorption places at the 02-terminated Fe304(lll) 
surface in a top and side view. 



in the most preferred site C is the same (-0.09/is) for 
GGA and GGA+J7 calculations. 



C. Adsorption on the O2 terminated surface 

For the oxygen termination the following four different 
adsorption sites (Fig. ^ were considered: the hollow over 
the Feoct2 atom (site A), a deep hollow over the Oi atom 
(site B), a hollow over the Fetot2 atom (site C), and the 
site above (on-top) the lower O2 (site D). All of the con- 
sidered positions appeared to be stable during structure 



optimization, and the energetically preferred site, both 
for Au and Pd adsorption, is site B, where the adsorbate 
atom is coordinated by three O atoms. 

The adsorbate induced changes in the relaxation pat- 
tern resulting from GGA and GGA+U calculations are 
qualitatively and quantitatively similar for both types of 
adsorbates (Fig. ITO)) . The adsorption of a Au/Pd atom 
on the O2 terminated surface suppresses the relaxations 
of the first two interplanar distances. The contraction of 
the first interplanar distance is reduced by one-half with 
respect to the clean surface, but the distance still remains 
25% shorter than that in the bulk crystal. The changes of 
the relaxations of deeper distances are rather small. In- 
terestingly, the difference in the sign of the relaxation of 
the fourth interlayer distance in the GGA and GGA+U 
results is removed in the adsorbate-oxide system. 

The adsorption binding energy in the considered sites 
and the adatom distance to the topmost oxide layer are 
displayed in Fig. [TlJ Additionally, Table IIIII shows nu- 
merical values for the most stable site. In general, the 
GGA+U binding is up to about 2 eV larger than that 
calculated within GGA. Similarly as for the Fctoti termi- 
nation the Pd atoms bind stronger than Au to the 02- 
terminated surface. The Au binding energy is 1.74 eV 
(GGA) and 3.66 eV (GGA+U), whereas for Pd it is 3.39 
eV (GGA) and 4.87 eV (GGA+C/). The bond lengths 
between Au and the three coordinating oxygens are not 
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FIG. 10: (Color online) Interlayer relaxation (in %) of the in- 
terplanar distance after Au and Pd adsorption in most stable 
sites on the O2 terminated Fe304(lll) surface. 
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FIG. 11: (Color online) The same as in Fig. [7] but for the O2 
terminated surface. 



symmetric and are 1.99, 2.07, and 2.07 A for GGA+U 
(2.03, 2.29 and 2.29 A for GGA). These nmnbers are very 
close to the corresponding Pd-0 bond lengths of 1.95, 
2.10, and 2.10 A. 

Figure [12] displays the layer-decomposed LDOS for the 
Au and Pd adsorbed in the energetically favored site B. 



TABLE 111: Adsorption energy, E^ds, adatom-surface dis- 
tance, dx-02, s-nd work function change, A"l>, for X (=Au,Pd) 
adatom in most stable sites on the 02-terminated surface. For 
each quantity the left and right hand side columns display 
GGA and GGA+U values, respectively. 



X (site) 


£ads (eV) 


dx-Oa 


(A) 


A* (eV) 


Au (B) 


1.74 3.67 


1.18 


0.39 


-1.91 -1.82 


Pd (B) 


3.39 4.87 


0.93 


0.82 


-2.05 -2.43 



The adsorbate covered surface remains half-metallic both 
when calculated in GGA and GGA+U although in the 
latter case a nonvanishing LDOS of the hybridized Au 
5d and O 2p states at the Fermi level makes the bands 
look more metallic-like. The Au and Pd states are much 
more delocalized and extended over a wider range of en- 
ergies compared with the Fctoti termination. The half- 
metallicity of the GGA-|-f7 bands of the Au/02 system is 
due to minority spin band of the Fe 3d states. The energy 
gap in the majority spin band is shifted to lower energies 
compared to the clean surface. Upon Au adsorption the 
hybridization of the Au 5d and O 2p states is strongest 
up to about 0.5 eV below the Fermi level while for the Pd 
4d states the hybridization with the O 2p states occurs 
up to the Fermi level. This stronger (compared to that of 
Au) hybridization of the oxide and Pd states at energies 
closer to the Fermi level explains the stronger binding 
of Pd than Au to the O2 terminated surface (Fig. [TT]) . 
The resulting LDOS are dominated by the O 2p electron 
states of the surface O2 layer hybridized with the Au 5d 
or Pd 4d states. The states of Fe atoms of the underlying 
Fetct2 and Fcocti layer contribute significantly only to the 
LDOS at the lower energy range, below —6 eV. For Pd 
adsorption the LDOS on the Pd atom is more localized 
than on Au and its main weight extends between —1.5 
eV and the Fermi level. In contrast to the clean termi- 
nation, in the presence of adsorbed Pd, responsible for 
half-metallic character of the surface is the energy gap in 
the minority spin band. 

The work function of the oxygen terminated surface de- 
creases dramatically upon Au/Pd adsorption (Fig. [TlT) . 
Au adsorption in the most preferred site B decreases 
the work function by about 1.9 eV (GGA) and 1.73 
eV {GGA+U). For Pd adsorption this decrease is even 
larger: 2.0 eV (GGA) and 2.4 eV (GGA+U) which means 
a nearly 30% reduction with respect to the clean surface 
value. A reduction of work function is connected with 
a decrease of the surface dipole moment and indicates 
a decrease in polarity. Both for Au and Pd a charge 
transfer from the adsorbate to the substrate is obtained. 
This is different from adsorption on the iron terminated 
surface, where a transfer from adsorbate to the surface 
is observed only for Pd, but is consistent with the large 
work function decrease described above. The Au adatom 
loses 0.70e (GGA) or 0.77e (GGA+U), Pd loses 0.89e 
(GGA) or 0.92e (GGA+U). Most of this charge is trans- 
ferred to the atoms of the O2 layer closest to the adatom 
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FIG. 12: (Color online) Same as in Fig. |8] but for Au/Pd adsorption on the O2 termination. 



(0.08e (GGA) and 0.30e (GGA-^t/) for Au, and 0.28e 
(GGA) and 0.25e (GGA+C/) for Pd adsorption). 

The magnetic moments on the oxygen atoms in the 
topmost surface layer are only little affected by the pres- 
ence of adsorbate. For a Au atom in the most pre- 
ferred site B, the GGA-I-C/ moments on the 0-atoms are 
changed by 0.11-0.16/LtB. A Pd atom in the same site 
induces the smaller changes (0.02-0.07/iB). For both ad- 
sorbates in the most stable site B, the moment of one of 
the O atoms in the topmost substrate layer is oriented 
down which is opposite to the direction of the remaining 
three O atoms, and is similar to the clean O2 terminated 
surface. The changes on the Fetet2 atoms of the first sub- 
surface layer resulting from the GGA-I-/7 calculations for 
Au and Pd adsorption are -0.40/xb and -0.29/x_b, respec- 
tively. A decreased magnitude of the magnetic moments 
of the Fetet2 atoms makes them comparable with those in 
the bulk magnetite. The changes of moments of atoms in 
deeper layers compared with the bulk are negligible. On 
the 02-terminated surface the magnetic moments on the 
Au differ from those on Pd atoms. The moments of the 
Au atom are lower than 0.1/is. Only the magnetic mo- 
ments of the Au atom in the site B and C resulting from 
the GGA are larger than that (O.lS/Lts and — 0.15^s, re- 
spectively). The moment on Au in the most preferred 
site B is — 0.02/^^ (GGA+C/). The magnetic moments 
on Pd are much larger than those on the Au and are 
in the range 0.4-0.9/iB. In the most stable site B the 



magnetic moment on Pd is O.Sl/is (GGA) and 0.92/iB 
(GGA+C/). 



IV. SUMMARY AND CONCLUSION 

We have presented a detailed DFT and DFT+C/ study 
of the structural, electronic, and magnetic properties of 
the clean magnetite (111) surface, and the adsorption 
of Au and Pd atoms on its two stable, Fctoti and O2, 
terminations. Inclusion of on-site Coulomb correlations 
in GGA+C/ approach modifies profoundly the electronic 
structure of magnetite surfaces. Based on ah initio ther- 
modynamics the Fctoti terminated surface is confirmed 
to be the most stable one over a broad range of oxygen 
pressures. It shows metallic character when calculated 
within GGA and half-metallic DOS with zeroth LDOS at 
the Fermi level when calculated within GGA+C/. All ter- 
minations studied exhibit a large inward relaxation of the 
first interlayer distance. Both Au and Pd bind strongly 
to the magnetite surface and induce large changes in the 
surface geometry. At the Fctcti termination different sites 
are favored for adsorption of Au and Pd. For Au adsorp- 
tion the most favorable site is on top of the Fctcti atom, 
whereas for Pd it is threefold coordinated hollow. At the 
O2 termination the threefold hollow site B, where the 
adatom is coordinated by three O atoms, is most stable 
both for Au and Pd adsorption. The GGA+C/ bonding 
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is respectively 0.2-0.5 eV and 0.5-1.5 eV stronger on the 
iron and the oxygen terminated surface, than that re- 
sulting from standard GGA calculations. The binding 
is stronger for Pd than Au and for both adsorbates is 
distinctly stronger on the oxygen than on the iron ter- 
minated surface. The Au/Pd bonding either to iron or 
oxygen terminated (111) surface of magnetite shows close 
resemblance to that reported by us for Au/Pd adsorption 
on the hematite (0001) surface terminations both with re- 
spect to the preference of the sites and strength of the 
bonding [32j |. 
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